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Abstract 

In this work we develop a tractable structural model with analytical default 
probabilities depending on a random default barrier and possibly random volatil- 
ity ideally associated with a scenario based underlying firm debt. We show how 
to calibrate this model using a chosen number of reference Credit Default Swap 
(CDS) market quotes. In general this model can be seen as a possible exten- 
sion of the time- varying ATIP model in Brigo and Tarenghi (2004). The calibra- 
tion capability of the Scenario Volatility /Barrier model (SVBATIP), when keeping 
time-constant volatility, appears inferior to the one of ATIP with time- varying de- 
terministic volatility. The SVBATIP model, however, maintains the benefits of 
time- homogeneity and can lead to satisfactory calibration results, as we show in a 
case study where we compare different choices on scenarios and parameters. 

Similarly to ATIP, SVBATIP is suited to pricing hybrid equity/credit deriva- 
tives and to evaluate counterparty risk in equity payoffs, and more generally to 
evaluate hybrid credit /equity payoffs. We consider the equity return swap in Brigo 
and Tarenghi (2004) and show its valuation under SVBATIP with the same CDS 
and equity calibration input used earlier for ATIP, and further we hint at equity 
default swap valuation in the conclusions. 



*We are grateful towards Eymen Errais, for suggesting us to have a look at the random barrier 
approach, and to Kay Giesecke for helpful discussion and correspondence. 
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Duke of Alencon, this was your default, 
That, being captain of the watch to-night. 
Did look no better to that weighty charge. 

King Henry VI, Part 1, Act 2, Scene I. 



1 Introduction 

Modelling firms default is becoming more and more important, especially in recent times 
where the market is experiencing a large development in credit derivatives trading. In 
this paper we develop a tractable structural model with analytical default probabilities 
depending on some dynamics parameters and on a random default barrier ideally as- 
sociated with the underlying firm debt. The scenario based default barrier models our 
possible uncertainty on the debt level of the firm. This uncertainty may be due to the 
fact that, for example, we do not trust completely the information from Balance Sheets. 

We show how to cahbrate the model using a chosen number of Credit Default Swap 
(CDS) market quotes. We apply the structural model to a concrete calibration case. 
This model can be seen as a variant or even a possible extension of the time-varying 
Analytically Tractable 1st Passage model (ATIP) in Brigo and Tarenghi (2004). There 
are a number of approaches for extending this earlier work. In Brigo and Tarenghi (2004) 
we essentially used information from equity volatility and CDS quotes to determine a 
time- varying value of the deterministic firm- volatility and a reference level for the default 
barrier implying exact calibration of a given number of CDS quotes. In this paper we 
consider a scenario based default barrier/volatility version of the basic model, which 
is a step forth, but take time-constant volatilities (in each scenario), which is a step 
backwards. Notwithstanding the time-constancy of the volatility scenarios, which is 
a restriction with respect to the earlier deterministic time-varying coefficients ATIP 
model, we refer to the random coefficients model as to an "extension", that we term 
Scenario Volatility/Barrier Analytically Tractable 1st Passage model (SVBATIP). The 
calibration of the CDS term structure in SVBATIP is left to the random default barrier, 
to the volatility scenarios parameters and to the probabilities of different scenarios. These 
can all be considered as calibrating parameters, or else some can be fixed arbitrarily and 
the cahbration can be left to the remaining ones. We illustrate different possible choices 
with a case study. 

With the first extension where we use scenarios only on the barrier (taking instead 
the same volatility in all scenarios) the model is not very flexible and can calibrate a 
low number of CDS, typically 3 or 4. Indeed, default barrier scenarios appear not to 
be natural fitting parameters for a large number of CDS quotes. In this case it can be 
best to associate scenarios also to volatilities, or to take a time varying volatility to be 
partially fitted to CDS quotes even in the case of random barrier. In this note we explore 
the former solution, and while time-constant volatility scenarios do not allow a perfect 
fit of all the five CDS maturities, they still allow a small cahbration error, which can be 
deemed to be acceptable when compared to the market induced bid ask spreads. What 
is more, having a model like SVBATIP with constant parameters, even if under different 
scenarios, may help with respect to robustness, typically endangered by all-fitting time- 
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varying functions as in the basic ATIP model. Nonetheless, the time-varying volatility 
scenarios extension will be addressed in further work. 

Finally we notice that, as for the earlier deterministic barrier ATIP model, the CDS 
cahbrated SVBATIP model is ideally suited to price hybrid equity/credit derivatives 
and to evaluate counterparty risk in equity payoffs. Given the same CDS and equity 
calibration inputs, we compare the price of an equity return swap under the deterministic 
time varying volatility model ATIP with the same price under the scenarios model 
SVBATIP addressed in this paper. 



2 Calibrating the deterministic barrier ATIP model 

The fundamental hypothesis of the model we resume here is that the underlying firm 
value process is a Geometric Brownian Motion (GBM), which is also the kind of process 
commonly used for equity stocks in the Black Scholes model. 

Classical structural models (Merton (1974), Black Cox (1976)) postulate a GBM 
(Black and Scholes) lognormal dynamics for the value of the firm V. This lognormal- 
ity assumption is considered to be acceptable. Crouhy et al (2000) report that "this 
assumption is quite robust and, according to KMVs own empirical studies, actual data 
conform quite well to this hypothesis.". 

In these models the value of the firm V is the sum of the firm equity value S and of 
the firm debt value D. The firm equity value S, in particular, can be seen as a kind of 
(vanilla or barrier-like) option on the value of the firm V. Merton typically assumes a 
zero-coupon debt at a terminal maturity T. Black Cox assume, besides a possible zero 
coupon debt, safety covenants forcing the firm to declare bankruptcy and pay back its 
debt with what is left as soon as the value of the firm itself goes below a "safety level" 
barrier. This is what introduces the need for barrier option technology in structural 
models for default. 

In Brigo and Tarenghi (2004) the following proposition is proved: 

Proposition 2.1. (Analytically- Tractable First Passage (ATIP) Model) Assume 
the risk neutral dynamics for the value of the firm V is characterized by a risk free rate 
r{t), a payout ratio q{t) and an instantaneous volatility cr(t), according to 

dV{t) = V{t) (r(t) - q{t)) dt + V{t) a{t) dW{t) 

and assume a safety barrier H{t) of the form 

H(t) = Hexp £ (^q(s) - r{s) + (1 + 2/3)^) ds^ (1) 

where j3 is a parameter that can be used to shape the safety barrier, H is a reference 
level, and let r be defined as the first time where V hits the safety covenants barrier H 
from above, starting from Vq> H, 

T = mi{t > : V(t) < H(t)}. 

Then the survival probability is given analytically by 



Q{t >T} = 
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In presence of zero coupon debt for the maturity T, the actual default time is r if 
T < T, and is T if r > T and the value of the firm at T is below the debt face value. 
Otherwise, if r > T and the value of the firm at T is above the debt face value, there is 
no default. 

In Brigo and Tarenghi (2004) we cahbrate CDS quotes by means of the above formula 
inserted in CDS valuation, backing out t i-^ a{t) from increasing-maturity CDS quotes 
through a cascade inversion method. We now recall briefly the CDS payoff and its risk 
neutral pricing formula. 

One of the most representative protection instruments that can be used against de- 
fault is the Credit Default Swap (CDS). Consider two companies "A" (the protection 
buyer) and "B" (the protection seller) who agree on the following. 

If a third reference company "C" {the reference credit) defaults at a time tc € (Ta, Tf,], 
"B" pays to "A" at time t = tc itself a certain "protection" cash amount Lgd (Loss 
Given the Default of "C"), supposed to be deterministic in the present paper. This 
cash amount is a protection for "A" in case "C" defaults. A typical stylized case occurs 
when "A" has bought a corporate bond issued from "C" and is waiting for the coupons 
and final notional payment from this bond: If "C" defaults before the corporate bond 
maturity, "A" does not receive such payments. "A" then goes to "B" and buys some 
protection against this risk, asking "B" a payment that roughly amounts to the bond 
notional in case "C" defaults. Typically Lgd is equal to a notional amount, or to a 
notional amount minus a recovery rate. We denote the recovery rate by "Reg" . 

In exchange for this protection, company "A" agrees to pay periodically to "B" a 
fixed "running" amount /?, at a set of times {Ta+i, . . . , T5}, = Ti — Ti_i, Tq = 0. These 
payments constitute the "premium leg" of the CDS (as opposed to the Lgd payment, 
which is termed the "protection leg"), and R is fixed in advance at time 0; the premium 
payments go on up to default time r if this occurs before maturity I],, or until maturity 
Th if no default occurs. 

"B" protection Lgd at default tc if Ta<rc <Tb "A" 

"B" ^ rate R at T„+i, . . . , T,, or until default tc ^ "A" 

Formally, we may write the RCDS ( "R" stands for running) discounted value at time 
t seen from "A" as 

b 

URCDSa,b{t) := -D{t,T){T - r^(^)_i)i?l{T,<r<r6} " ^ D{t,Ti)aiRl{r>T,} + 

i=a+l 

'i-{Ta<T<n}D{t,T)LGJ3 (3) 

where t e [T^(f)_i, T^(()), i.e. T^(t) is the first date among the Tj's that follows t, and where 
ctj is the year fraction between T'j_i and Tj. The pricing formula for this payoff depends 
on the assumptions on the interest rates dynamics and on the default time r. Let J-'t 
denote the basic filtration without default, typically representing the information fiow of 
interest rates and possibly other default-free market quantities (and also intensities in 
the case of reduced form models), and — ^tVcx ({r < u},u < t) the extended filtration 
including explicit default information. In our earlier "structural model" ATIP framework 
with deterministic default barrier the two sigma-algebras coincide by construction, i.e. 
Qt — ^t, because here the default is completely driven by default-free market information. 
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This is not the case with intensity models, where the default is governed by an external 
random variable and is strictly included in Qt, i.e. J-'t C Qf 

If we include barrier scenarios and take the barrier independent of W, we are back in 
a situation where J-'t <Z Qf Indeed, if for example we consider the random barrier version 
with deterministic volatility and we know that default occurs at a given time, say now, 
we also know that V has hit one of the barriers in the current scenario. By observing V, 
we know which barrier scenario has realized itself, and by knowing that default is now 
we know that the value of V now is telling us the value of H. But observation of V alone 
(as happens with J^) until "now" would not have allowed us to know H now or to know 
that default is "now" . Thus, contrary to the basic ATIP model, Qt tells us more than J^t 
under scenario based default barrier. The situation becomes more involved if we allow 
for a random volatility a, independent of all observable quantities. We might assume the 
random volatility is realized after a small time e, as we do in smile modeling for example 
in Brigo, Mercurio and Rapisarda (2004), but this leads to undesirable features such as 
unrealistic future dynamics. While in smile modeling keeping a random "forever" would 
have led us to an underlying equity or FX rate that is not observed, unacceptably leading 
to payoffs that remain random even at maturity, in case of structural models we may 
argue that the current value of the firm is not necessarily known with certainty. This 
would cut us some slack in keeping the volatility in V random "forever" . The mechanics 
of iterated expectation in computing prices does not change, one just has to make sure to 
increase the inner conditioning filtration with the unobservable variable and everything 
works fine. Here, however, even knowledge of r is not necessarily implying knowledge 
of H or V, given V^s uncertainty. At r wc would only know that V hit H, but since 
V is now random, depending on the volatility scenario, we cannot be sure of which H 
scenario has been hit. Thus the Q filtration is less informative in this random volatility 
case. 

More generally, a discussion on the role of filtrations under random default barriers 
and partial observation is in many papers by Giesecke, see for example Giesecke and 
Goldberg (2004) (who introduce the P model). 

We denote by CDSa,b(i, R, Lgd) the price at time t of the above standard running 
CDS. In general we can compute the CDS price according to risk-neutral valuation (see 
for example Bielecki and Rutkowski (2001)): 

CDS,,b(t,i?,LGD) = E{Unc^Sa,b{t)\gt} =: Et{Uncusa,b{t)} (4) 

in our structural model setup. Straightforward computations lead to the price at initial 
time as 

CDSa,6(0, R, Lgd) =R [ " PiO, t) {t - Tp^t)-i)dQiT > t) (5) 
-R ^ P(0, T,)aiQ{T > Ti) - Lgd /" ' P(0, t)dQ(r > t) 

i=a+\ "^^^ 

SO that if one has a formula for the curve of survival probabilities t h- > Q(r > t), as in 
our ATIP structural model, one also has a formula for CDS. 

A CDS is quoted through its "fair" R, in that the rate R that is quoted by the market 
at time t satisfies CDSa,6(t, R, Lgd) = 0. The fair rate R strongly depends on the default 
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probabilities. The idea is to use quoted values of these fair i?'s with different maturities 
to derive the default probabilities assessed by the market. 

While in simple intensity models the survival probabilities can be interpreted as 
discount factors (with credit spreads as discounting rates), and as such can be easily 
stripped from CDS's or corporate bonds, in structural models the situation is much 
more complicated. In fact, here, it is not possible to find a simple "credit spread" 
formulation for dQ{T > t} starting from ([2]). 

To calibrate the ATIP model to quoted market i?'s for different CDS, we insert the 
quoted i?'s in Formula ([5]) and find the a and H values that set said formula to zero. See 
Brigo and Tarenghi (2004) for several numerical examples and a case study on Parmalat 
CDS data. 

2.1 Numerical example 

In this section we present some results of the calibration performed with the ATIP struc- 
tural model. We borrow from Brigo and Tarenghi (2004). We consider CDS contracts 
having Vodafone as underlying with recovery rate Rec = 40% (Lgd = 0.6). In Table [T] 
we report the maturities Tj, of the contracts and the corresponding "mid" CDS rates 
i?Q™(0) (quarterly paid) on the date of March 10th, 2004, in basis points {Ibp = 10~^). 
We take = in all cases. 





CDS maturity 


i?o7(0) (bps) 


Kim 




ly 


20-mar-05 


19 


24 


21.5 


3y 


20-mar-07 


32 


34 


33 


5y 


20-mar-09 


42 


44 


43 


7y 


20-mar-ll 


45 


53 


49 


lOy 


20-mar-14 


56 


66 


61 



Table 1: Maturities of quoted CDS's with their corresponding spreads on March 10, 2004. 



CDS mat 


CDS value bid (bps) 


CDS value ask (bps) 


ly 


2.56 


-2.56 


3y 


2.93 


-2.93 


5y 


4.67 


-4.67 


7y 


24.94 


-24.94 


lOy 


41.14 


-41.14 



Table 2: CDS values computed with deterministic default intensities stripped from mid R 
quotes but with bid and ask rates R in the premium legs. 

In Table [2] we report the values (in basis points) of the CDS's computed inserting 
the bid and ask premium rate R quotes into the payoff and valuing the CDS with 
deterministic intensities stripped by mid quotes. This way we transfer the bid offer 
spread in the rates i? on a bid offer spread on the CDS present value. In Table [3] we 
present the results of the calibration performed with the structural model and, as a 
comparison, of the calibration performed with a deterministic intensity (credit spread) 
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model (using piecewise linear intensity). In this first example the parameters used for 
the structural model have been selected on qualitative considerations, and are /? = 0.5 
and H/Vq = 0.4 (this is a significant choice since this value is in line with the expected 
value of the random H, completely determined by market quotes, in the Scenario based 
model presented later on). 

We report the values of the calibrated parameters in the two models (volatilities and 
intensities) and the survival probabilities. 



n 


a 


Surv. 


Intensity 


Surv. 





36.625% 


100.000% 


0.357% 


100.000% 


ly 


36.625% 


99.627% 


0.357% 


99.627% 


3y 


17.311% 


98.316% 


0.952% 


98.316% 


5y 


17.683% 


96.355% 


1.033% 


96.355% 


7y 


17.763% 


94.206% 


1.189% 


94.206% 


lOy 


21.861% 


89.650% 


2.101% 


89.601% 



Table 3: Results of the cahbrations performed with both models. 

Comments on the realism of short term credit spreads and on the robustness of default 
probabilities with respect to CDS are in Brigo and Tarenghi (2004). 

The parameter P has not a direct economic meaning and in principle can be chosen 
arbitrarily, and can help in setting more or less stringent safety covenants, depending on 
one's attitude towards default risk. As for H, in Brigo and Tarenghi (2004) we present 
some methods to estimate it based on market considerations. 

3 The scenario volatility /barrier SVBATIP model 

In some cases it can be interesting to keep the volatility of the process V as an exogenous 
input coming from the equity and debt worlds (for example it could be related to an 
historical or implied volatility). Or we might retain a time- varying volatility to be used 
only partly as a fitting parameter. Or, also, we might wish to remove time-varying 
volatility to avoid an all-fitting approach that is dangerous for robustness. 

In all cases we would need to introduce other fitting parameters into the model. One 
such possibility comes from introducing a random default barrier. This corresponds to 
the intuition that the balance sheet information is not certain, possibly because the com- 
pany is hiding some information. In this sense, assuming a random default barrier can 
model this uncertainty. This means that we retain the same model as before, but now 
the default barrier level H is replaced by a random variable H assuming different scenar- 
ios with given risk neutral probabilities. At a second stage, we may introduce volatility 
scenarios as well. By imposing time-constant volatility scenarios, in each single scenario 
we loose ficxibility with respect to ATIP, but regain fiexibility thanks to the multiple 
scenarios on otherwise too simple time-constant volatilities. Even so, the scenario based 
model results in a less fiexible structure than the old deterministic volatility and barrier 
ATIP model as far as CDS cahbration is concerned. 
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In detail: 

Definition 3.1. (Scenario Volatility and Barrier Analytically Tractable 1st 
Passage model, SVBATIP) Let the firm value process risk neutral dynamics be given 
by 

dV{t) = (r(t) - q{t))V{t)dt + u{t)V{t)dW{t), 

same notation as earlier in the paper. This time, however, let the safety barrier param- 
eter H in (Cp and the firm value volatility function t ^— i/{t) assume scenarios 
{Hi,t I— > cr^(t)), . . . , {Hn-i, t a^~^(t)), {Hn, t I— > cr^(t)) with Q probability pi, . . . ,pn-i,Pn 
respectively. The safety barrier will thus he random and equal to 

H\t) := exp £ [q{s) - r{s) + (1 + 2/5)^^^^ ds^ 

with probability Pi. The random variables H and v are assumed to be independent of the 
driving Brownian motion W of the value of the firm V risk neutral dynamics. 

This definition has very general consequences. Indeed, if we are to price a payoff 11 
based on V, by iterated expectation we have 

N 

E[n] = E{E[U\H, z/]} = J2pM^\H = Hi,iy = a'] 

i=l 

Now, thanks to independence, the term E[n|if = Hi, v = a*] is simply the price of the 
payoff n under the model with deterministic barrier and volatility seen earlier in the 
paper, when the barrier parameter H is set to Hi and the volatility to t i— *• cr*(t), so that 
the safety barrier is H^. This means that, in particular, for CDS payoffs we obtain 

N 

CDSa,fe(t, R, Lgd) = J2 CDS„,b(t, R, Lgd; Hi, a') ■ p, (6) 

where CDSa,f,(t, R, Lgd; Hi, cr*) is the CDS price ([5]) computed according to survival prob- 
abilities ([2]) when the barrier H is set to Hi and the volatility to a*. Let us consider now 
a set of natural maturities for CDS quotes. This is to say that we assume = and T;, 
ranging a set of standard maturities, T;, = ly, 3y, 5y, 7y, lOy. Let us set 

CDS^ := CDSo,fe(0,i?,LGD;/f„(T,), 

i.e. the CDS with first reset in 0, final maturity and valued under the deterministic 
barrier model with the barrier set to if* and the deterministic volatility set to t \—>- cr*(t), 
which from now on we assume to be constant = a* for all t. 

Now assume we aim at calibrating the scenario based barrier/volatility model to a 
term structure of CDS data. Our first attempt is only with a scenario based barrier 
model, in that we take a common, time constant volatility a* = (x in all scenarios, so 
that the only uncertain quantity in the model is the safety level barrier parameter H. 
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3.1 CDS Calibration: Scenario based Barrier with Linear Al- 
gebra 

For this particular case we could reason with some linear algebra. Assume we have N 
CDS maturities to calibrate and consider the parameters we might use. If we resort 
to a scenario based barrier and fix the volatility exogenously, with barrier scenarios 
we have N values for H and — 1 values for the probabilities p, with a total of 2N — 1 
parameters to fit A^ quotes. Thus we have a system with more unknowns than equations. 
One possible way of proceeding would be to fix A^ — 1 parameters (some barriers and 
some probabilities) and then looking for the remaining free parameters. This procedure, 
when tried, is not very stable and the calibration is quite difficult. A better way to face 
the problem is the following. Let Ck^i ■= CDS^ = CDSo,fc(0, R, Lgd; Hi, a) 

If we write the relationship equating in correspondence of market i?'s to zero, as 
should be, for A^ maturities and with A^ barriers, we have a linear system of the form 

C-p = (7) 

where p is the (A^ x 1) column vector of the probabilities and is a (A^ x 1) null vector. It 
is well known from basic algebra that if the determinant of C is different from zero, then 
the system ([7]) admits a unique solution and this solution is identically equal to zero (the 
kernel of C is the trivial null space), which is not good as a set of probabilities p. Since 
we aim at a meaningful result for p we ask that det(C) = 0, so that the infinite solutions 
of the system will be given by all the eigenvectors relative to the null eigenvalue. Of 
all eigenvectors, we consider unitary norm versions. Hence the first aim is to find those 
values of the barriers that make the determinant of C equal to zero. This is a single 
equation with A^ unknowns. So we can fix A^ — 1 barrier parameters H and look for the 
last one solving this "vanishing determinant" equation. Obviously we must be careful 
in the choice of the parameters, since in principle given the first A^ — 1 barriers, the 
equation could admit no solutions. Also, even if we find an admissible value for the last 
barrier, there is no guarantee that the probabilities be all positive, since the eigenvectors 
in principle could have negative components. 

Actually, as one could guess beforehand, these considerations result in a calibration 
technique that does not seem to be very powerful. Let us analyze the particular case 
of the Vodafone telecommunication company at the date of March 10th, 2004 (the data 
are the mid quotes reported in Table [H the recovery rate is Reg = 40%, the common 
and time-constant volatility in all scenarios is a = 24%, (3 is set equal to 0.5 and the 
spreads are expressed in basis points). In the case of deterministic barrier we find a term 
structure of deterministic firm value volatility to calibrate exactly the data, as we have 
seen in detail in Section 12.11 

Here instead we keep the volatility fixed as an exogenous constant and look for the 
barrier parameters H and the probabilities p. 

First of all we try to calibrate the first two maturities. We have to fix the value for 
one of the barriers. We choose a high value, i.e. Hi = 0.8. The results of the calibration 
are reported in Table |H 

We see that the calibration assigns a low probability to the high barrier and a high 
probability to the low barrier: This is consistent with the fact that the CDS rates are 
very low, indicating a quite good credit quality of the underlying. In fact the high barrier 
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i 


Hi 


Pi 


1 
2 


0.3710* 
0.8000 


98.86% 
1.14% 



Table 4: Results of the calibration with the first two maturities. The asterisk indicates 
that the corresponding value is the value obtained from the calibration while the other 
barrier value has been fixed in advance. 

(corresponding to high debt) indicates proximity to default, while the low barrier (low 
debt) indicates a distant default possibility. 

Now let us add a maturity, i.e. let us move from two maturities to three maturities, 
namely the first three dates in Table [H It would be somehow natural to try a configura- 
tion of barriers which is an extension of the previous one, i.e. to fix the two barriers just 
found and to look for a third one. However, this procedure does not lead to a solution: 
To calibrate we need to change one of the barriers. We decided to keep fixed H = 0.8 
and change the other one, choosing H = 0.2. The results are reported in Table [51 



i 


H, 


Pi 


1 


0.2000 


61.06% 


2 


0.4303* 


37.83% 


3 


0.8000 


1.10% 



Table 5: Results of the calibration with the first three maturities. The asterisk indicates 
that the corresponding value is the value obtained from the calibration while the other 
values have been fixed in advance. 

We notice that a splitting has occurred: From a single barrier set at 0.3710 we 
obtained two barriers, one higher and one lower. The highest barrier has still a low 
probability, almost equal to the corresponding one obtained earlier, when calibrating 
only two maturities. The other probability of the two- maturities case has split, following 
the splitting of the H parameters, and now the lowest H is the most likely. 

Again, moving to the first four maturities in Table [H we would try simply to add a 
new barrier parameter while retaining the if's calibrated under three quotes. As before, 
this calibration does not work. Then we substitute the last barrier found (0.4303) with 
0.5500 before trying the calibration. In Table [6] we report the result of the calibration. 



i 


Hi 


Pi 


1 


0.2000 


27.65% 


2 


0.3347* 


66.75% 


3 


0.5500 


7.71% 


3 


0.8000 


0.89% 



Table 6: Results of the calibration with the first four maturities. The asterisk indicates 
that the corresponding value is the value obtained from the calibration while the other 
values have been fixed in advance. 

As when moving from two to three maturities, when moving from three to four 
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maturities we observe a splitting of one barrier parameter H into two values, one higher 
and one lower. We see that the greatest part of the probability is divided among the three 
"low" barriers while the highest barrier has still a low probability. More precisely, we see 
that adding information (i.e. adding maturities), the high barrier seems to become less 
probable. In general adding information leads obviously to a refinement of the barrier 
configuration. 

It may be interesting to see what happens in case of a worse credit quality. To see 
this let us arbitrarily change the CDS rates (all other things being equal), in particular 
let us double them, as shown in Tabled and repeat the calibration as before. The results 
are reported in Tables M and M 



Maturity Tf^ 


Doubled Rate i?^|^°(0) (bps) 


March 21st, 2005 
March 20th, 2007 
March 20th, 2009 


43.0 
66.0 
86.0 



Table 7: Maturities of quoted CDS's with their corresponding doubled spreads, valuation 
date March 10, 2004. 



i 


Hi 


Pi 


1 

2 


0.4123* 
0.8000 


97.76% 
2.24% 



Table 8: Results of the calibration with the first two maturities. The asterisk indicates 
that the corresponding value is the value obtained from the calibration while the other 
barrier value has been fixed in advance. 



i 


Hi 


Pi 


1 


0.2000 


19.81% 


2 


0.4277* 


77.98% 


3 


0.8000 


2.21% 



Table 9: Results of the calibration with the first three maturities. The asterisk indicates 
that the corresponding value is the value obtained from the calibration while the other 
values have been fixed in advance. 

We see that by increasing the CDS rates, i.e. decreasing the credit quality of the 
underlying firm, the calibrated barriers parameters H become higher and also the prob- 
abilities to have high barriers increase, as one would expect, since high barriers increase 
the proximity of default. 

With this method we have not been able, by adding a further scenario to H, to 
calibrate all of the five CDS of Table [H This is due probably to the time-homogeneity of 
the model. To have a larger calibration capability we may have to use a time varying u. 

To sum up, as one could have expected, this "linear algebra" procedure and the 
scenario based barrier approach in general are not very robust and are quite sensitive to 
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CDS rates: For high rates the cahbration works only for a low number of maturities, and 
the refinement is limited. This is one of the main drawbacks of the model which makes 
it less powerful than the standard calibration with deterministic barrier ATIP of Brigo 
and Tarenghi (2004), at least if one insists in having an exact calibration. One possibility 
to improve the model is to include a time dependency for the parameters (for example 
in the volatility a), as hinted at before. This idea is currently under investigation. A 
different possibility is moving in two different directions: first, abandoning the "linear 
algebra" approach and using a numerical optimization within the same framework, and 
secondly abandoning the single time-constant volatility scenario and allowing several 
time-constant volatility scenarios to come into play, which we explore in Section 13.31 

3.2 CDS Calibration: Scenario based Barrier with Numerical 
Optimization 

Here we abandon the requirement to obtain a perfect calibration and allow for a (small) 
calibration error, while trying to restore a one to one correspondence between model 
parameters and market quotes. 

We start by the first three quotes in Table [H but this time calibrate them with 
an optimization, abandoning the "linear algebra" approach. We minimize the function 
expressing the sum of the squares of the CDS's prices in the model, since each CDS 
price should be zero in correspondence of the market quoted R. Thus we actually solve 
numerically 

3 

[HIH;,pI] = argmin^,p^[piCDSo,fc(0,i?oT(0),LGD;i/i) 

fc=i 

+ (1 - pi)CDSo,fc(0, i?oT(0)' Lgd; 

We find the following results: 





Pi 


0.3188 
0.6592 


94.83% 
5.17% 



Table 10: SBATIP calibrated parameters. Exact calibration of the first three CDS 
quotes of Table [H The expected value of the barrier is E[if] = 0.3364. 

Luckily, in this case the target function is practically zero, so that the calibration is 
exact. While this could be somehow hoped, since the number of free parameters matches 
the number of quotes to be calibrated, in general it is not guaranteed. Indeed, we will 
see shortly that even if the parameters are in the right number they often do not have 
enough flexibility to account for the market quotes exactly. 

Let us move to optimizing over the five quotes of Table [U To this end, we can choose 
the version with three H scenarios and two probabilities. We run an optimization to 
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solve 

5 



k=l 



PiCDSo,fc(0,i?S;^°(0),LGD;i7i 



where in the optimization we impose the p to take values that are allowed for probabili- 
ties, i.e. in [0, 1] and adding up to one. This time we obtain 





Pi 


0.7296 
0.3384 
0.7296 


1.24% 
97.52% 
1.24% 



Table 11: SBATIP calibrated parameters. Calibration of the first five CDS quotes of 
Table [H The expected value of the barrier is K[H] = 0.3481. 



The calibration error is not zero this time, but amounts to 9156ps^. We will discuss 
the calibration error shortly. One barrier scenario H remains close to 0.3188, with an even 
higher probability than before. Notice an important point: the two other scenarios of 
H more in proximity of default are identical, being both 0.7296, with equal probabilities 
0.0124. This suggests the parametrization to be not effective in describing a dynamics 
consistent with the cross sectional data we are observing, since the parameters collapse to 
a sub-parameterization with just two barrier scenarios (three parameters in total) when 
fitting five quotes. This seems to suggest the new parameters we added with respect to 
the previous case add nothing to explain the increase of data. And indeed, if one tries 
the model with two barriers and one probability to fit the five quotes above, consistently 
with what we just obtained one finds the same calibration error and the parameters 



barrier 


probability 


0.7296 
0.3384 


2.48% 
97.52% 



Table 12: SBATIP calibrated parameters. Calibration of the first five CDS quotes of 
Table [1] with only three unknowns parameters. 



Let us have a look at the calibration error in single CDS's. If we compute the five 
CDS prices using the R mid quotes and with the parameters H and p above coming out 
of the optimization, we obtain the values reported in Table [131 

We can compare these calibration errors on the single CDS's present values with the 
errors induced in the CDS present value by the market R bid-ask spreads, reported in 
Table [2l The first two CDS are out of the bid ask windows with our calibration, the 
worst case being the second, giving 9.99 bps of present value against a corresponding bid 
ask window [-2.93,2.93]. 

This situation cannot be considered to be fully satisfactory, which leads us to the 
following extension. 
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CDS maturity 


nF>Q {r\ DMID TT TT II . ^ ^ ^ \ {U-r^i-A 


ly 


-Z. 1 1 




Q QQ 

y .yy 


5y 


-1.47 


7y 


-22.99 


lOy 


16.63 



Table 13: CDS values obtained using the parameters resulting from the calibration. 



3.3 Scenario based Barrier and Volatility: Numerical Opti- 
mization 

Now we consider again all of the five quotes of Table [1] and we use the general model 
with two scenarios on the barrier/volatility parameters (H,!/) and one probability (the 
other one being determined by normalization to one), with a total of five parameters for 
five quotes. Here, by trial and error we decided to set /3 = 0. 



[Hl,H;;al,ai;pl 



argmmj:^^^ 



E 

fc=i 



PiCDSo,fc(0,i?ST(0),LGD;i7i,ai) 



-f(l - pi)CDSo,fc(0, K'l^iO), Lgd; H2, a' 





0"i 


Pi 


0.3721 
0.6353 


17.37% 
23.34% 


93.87% 
6.13% 



Table 14: SVBATIP model calibrated to Vodafone CDS. The expected value of the 
barrier is E[i7] = 0.3882. 



We obtain a much lower optimization error than before, i.e. 147&ps^, corresponding 
to the following calibration errors on single CDS present values: 



CDS maturity 


CDSo,fc(0, R^f, H,, H,; a\ a'; p,) (bps) 


ly 


1.38 


3y 


-3.89 


5y 


8.16 


7y 


-7.56 


lOy 


2.41 



Table 15: CDS values obtained using the parameters resulting from the calibration. 

Now the single CDS calibration errors are much lower than before, being the CDS 
present values corresponding to market R closer to zero, with the exception of the five 
years maturity, which could be adjusted by introducing weights in the target function. 
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In particular we introduce weights that are inversely proportional to the bid-ask 
spread and repeat the calibration for the SVBATIP model finding the results presented 
in Tables [16] and [T7| (compare with Tables [H] and [T5|) . 







Pi 


0.3713 
0.6239 


17.22% 
22.17% 


92.63% 
7.37% 



Table 16: SVBATIP model calibrated to Vodafone CDS using weights in the objective 
function which are inversely proportional to the bid ask-spread. The expected value of 
the barrier is E[H] = 0.3899. 



CDS maturity 


CDSo,fc(0, R^^°, Hu H2; a\ a^- p,) (bps) 


ly 


5.85 


3y 


-3.76 


5y 


4.92 


7y 


-10.46 


lOy 


1.47 



Table 17: CDS values obtained using the parameters resulting from the weighted cali- 
bration. 

We see that now the ly CDS value is outside the bid-ask spread, but conversely the 
5y CDS value has decreased, assuming a size more in line with the spread. This is what 
we aimed at, since the 5y CDS is probably the most liquid one and the model needs to 
reproduce it well. 



0.200% 7 



0.100% 




-0.200% ^ 



Figure 1: SBATIP-ATIP survival probabilities difference. The data for the SBATIP are those 
in Table [10] and refer to an exact calibration of 3 CDS quotes. 

Furthermore, to see the goodness of the calibration, we can compare the survival prob- 
abilities of the S(V)BAT1P model with those of the perfectly calibrated ATIP model. 
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0.400% 7 



0.200% 




-0.400% ^ 



Figure 2: SVBATIP-ATIP survival probabilities difference. The data for the SVBATIP are 
those in Table [16] and refer to an optimal minimum calibration of 5 CDS quotes using weights 
in the objective function. 

In Figure [T] we plot the difference between the survival probabilities computed with the 
SBATIP model exact calibration (see Table [TU]) and the ATIP calibration (on a 5 years 
horizon). In Figure [2] we plot the difference between the survival probabilities computed 
with the SVBATIP model weighted calibration (see Table [T6|) and the ATIP calibration 
(on a 10 years horizon). The two figures show a good agreement between the models in 
terms of probabilities, and the main difference is that the survival probabilities nearly 
coincide at CDS maturities in the case of exact calibration (Figure [1]), while a differ- 
ence is found in the minimization case (Figure [2]). Anyway, it is clear that the fit is 
quite good, indicating two main results: (i) the SVBATIP model returns a good fit of 
CDS quotes, even if the calibration is not perfect; (ii) as already observed in Brigo and 
Tarenghi (2004), CDS are a robust source for survival (default) probabilities, whatever 
model is used. 

More in general, when combining default barrier and volatility scenarios, scenarios on 
u and H can be taken jointly (as we did so far) or separately, but one needs to keep the 
combinatorial explosion under control. The pricing formula remains easy, giving linear 
combination of formulas in each basic scenario. Taking time-varying parametric forms 
for the (t"s can add flexibility and increase the calibrating power of the model, and will 
be addressed in future work. 

A final important remark is in order. We see that the parameters resulting from 
the scenario versions calibration are more credible than those obtained in the ATIP 
framework. In particular we notice that in all cases we have an expected H which 
is comparable with the fixed H = 0.4 used in the deterministic case. What is more, 
even if the H are similar, the volatilities involved are quite smaller than before. As 
previously hinted at, this fact is essentially due to volatility scenarios inducing a mixture 
of lognormal distributions for the firm value, implying fatter tails, and allowing for the 
same default probabilities with smaller volatilities. 
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4 Counterparty risk in equity swap pricing: Reprise 

"Of course it does. 

A conservative estimate puts material reality as ten thousand million light years across. 
Even using the combined computational function of all the processors in the united plan- 
ets capable of dealing with numerical scale, it would take at least ninety years to..." 

Brainiac 5, "Legion Lost", 2000, DC Comics 

This section summarizes the results on counterparty risk pricing in Equity Return 
Swaps under AT IP in Brigo and Tarenghi (2004) and presents a short analysis based on 
the SB VAT IP models introduced in this paper. This is an example of counterparty risk 
pricing with the calibrated structural model. 

Let us consider an equity return swap payoff. Assume we are a company "A" entering 
a contract with company "B", our counterparty. The reference underlying equity is 
company "C". The contract, in its prototypical form, is built as follows. Companies 
"A" and "B" agree on a certain amount K of stocks of a reference entity "C" (with price 

5 — 5"") to be taken as nominal (A^ = K Sq). The contract starts in = and has 
final maturity Ti, — T. At t — there is no exchange of cash (alternatively, we can think 
that "B" delivers to "A" an amount K of "C" stock and receives a cash amount equal to 
KSq). At intermediate times "A" pays to "B" the dividend flows of the stocks (if any) 
in exchange for a periodic rate (for example a semi-annual LIBOR or EURIBOR rate 
L) plus a spread X. At final maturity T = Tj,, "A" pays KSt to "B" (or gives back the 
amount K of stocks) and receives a payment KSq. This can be summarized as follows: 

Initial Time 0: no flows, or 
A — ^ KS^ cash — > B 
A< — K equity of "C" < — B 

Time Tf. 

A — > equity dividends of "C" — > B 
A i — Libor + Spread < — B 

Final Time T^: 
A — > K equity of "C" — > B 
A < — KS^ cash < — B 

The price of this product can be derived using risk neutral valuation, and the (fair) 

spread is chosen in order to obtain a contract whose value at inception is zero. We ignore 
default of the underlying "C" , thus assuming it has a much stronger credit quality than 
the counterparty "B". It can be proved that if we do not consider default risk for 
"B" either, the fair spread is identically equal to zero. But when taking into account 
counterparty default risk in the valuation the fair spread is no longer zero. In case an 
early default of the counterparty "B" occurs, the following happens. Let us call t = tb 
the default instant. Before r everything is as before, but if r < T, the net present value 
(NPV) of the position at time r is computed. If this NPV is negative for us, i.e. for "A" , 
then its opposite is completely paid to "B" by us at time r itself. On the contrary, if it 
is positive for "A" , it is not received completely but only a recovery fraction Rec of that 
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NPV is received by us. It is clear that to us ("A") the counterparty risk is a problem 
when the NPV is large and positive, since in case "B" defaults we receive only a fraction 
of it. 

The risk neutral expectation of the discounted payoff is given in the following propo- 
sition (see e.g. Brigo and Tarenghi (2004), L{S,T) is the simply compounded rate at 
time S for maturity T): 

Proposition 4.1. (Equity Return Swap price under Counterparty Risk). The 

fair price of the Equity Return Swap defined above can be simplified as follows: 

ERS{0) = KSoXJ2(^^P{^,Ti) - LanEollir<n}D{0,r){NPV{T))A. 

1=1 ^ J 

where 

NPV{t) = eJ-KNPVZJ^^^,,{t)+KSo J2 D{r,T,)adm-i,T.) + X) 

^ i=/3(r) 

+ {KSo-KST,)D{T,n)y (8) 

and where we denote by NPYi^^^^^^^^iu) the net present value of the dividend flows be- 
tween s and t computed in u. 

The first term in Hes is the equity swap price in a default-free world, whereas the 
second one is the optional price component due to counterparty risk. 

If we try and find the above price by computing the expectation through a Monte 
Carlo simulation, we have to simulate both the behavior of St for the equity "C" under- 
lying the swap, and the default of the counterparty "B" . In particular we need to know 
exactly r = r^. Obviously the correlation between "B" and "C" could have a relevant 
impact on the contract value. Here the structural model can be helpful: Suppose to cal- 
ibrate the underlying process V to CDS's for name "B", finding the appropriate default 
barrier and volatilities according to the procedure outlined earlier in this paper with the 
ATIP model. We could set a correlation between the processes V^^ for "B" and St for 
"C" , derived for example through historical estimation directly based on equity returns, 
and simulate the joint evolution of \Vt^ , St]. As a proxy of the correlation between these 
two quantities we may consider the correlation between S^ and S'f , i.e. between equities. 

Going back to our equity swap, now it is possible to run the Monte Carlo simulation, 
looking for the spread X that makes the contract fair. 

We performed some simulations under different assumptions on the correlation be- 
tween "B" and "C". We considered five cases: p = —1, p = —0.2, p = 0, p = 0.5 and 
p = 1. In Table [TS] we present the results of the simulation, together with the error 
given by one standard deviation (Monte Carlo standard error). For counterparty "B" 
we used the Vodafone CDS rates seen earlier. For the reference stock "C" we used a 
hypothetical stock with initial price 5*0 = 20, volatility a = 20% and constant dividend 
yield q = 0.80%. The contract has maturity T = 5y and the settlement of the LIBOR 
rate has a semi-annual frequency. Finally, we included a recovery rate Rec = 40%. The 
starting date is the same we used for the calibration, i.e. March 10th, 2004. Since the 
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reference number of stocks K is just a constant multiplying the whole payoff, without 
losing generality we set it equal to one. 

In order to reduce the errors of the simulations, we have adopted a variance reduc- 
tion technique using the default indicator (whose expected value is the known default 
probability) as a control variate. In particular we have used the default indicator l{r<r} 
at the maturity T of the contract, which has a large correlation with the final payoff. 
Even so, a large number of scenarios is needed to obtain errors with a lower order of 
magnitude than X. In our simulations we have used N = 2000000. 

We notice that X increases together with p, and in Brigo and Tarenghi (2004) we 
explain why this is natural. 



p 


X 


ES payoff 


MC error 


-1 











-0.2 


2.45 


-0.02 


1.71 





4.87 


-0.90 


2.32 


0.5 


14.2 


-0.53 


2.71 


1 


24.4 


-0.34 


0.72 



Table 18: Spread X (in bps) under five correlations, 5*0 = 20, basic ATIP model. We also re- 
port the value of the average of the simulated payoff (times 10000) across the 2000000 scenarios 
and its standard error, thus showing that X is fair (leads to a zero NPV). 

To check the impact of the scenarios barrier, we have re-priced with the same X's 
found in Table [18] for ATIP our equity swap under the SBATIP model calibrated to 
the same CDS data, i.e. with the parameters given in Table [THl We have the outputs 
given in Table [191 If the models are close in terms of equity swap pricing, this payoff 
should be zero. We find actually values ranging from —28bps to 165 bps. Recalling that 
S'o = 20, for a unit notional we would get a maximum discrepancy of about 165/20 ~ 8 
bps. By comparing with the bid-offer window induced by market bid-offer i?'s on CDS's 
with unit notional, as shown for example in Table [21 we realize that this difference, if not 
completely neghgible, is small. Besides, for the largest payoff value, i.e. p = 0.5 (where 
the ATIP model gives -0.53 bps and the SBATIP gives 165.5 bps) we notice that with 
X = 16 (instead of 14.2) the ATIP model would give a payoff price of 166.9, so that our 
price difference between SBATIP and ATIP, when translated back in ATIP spread X, 
is less than two basis points. 



p 


X 


ES payoff 


MC error estimate 


-0.2 


2.45 


-28.44 


1.49 





4.87 


3.45 


2.04 


0.5 


14.2 


165.50 


2.29 



Table 19: Equity swap valuation under the SBATIP model calibrated to the same CDS 
data as the ATIP model leading to Table [IS 

We can try the same considerations with the SVBATIP model, with the weighted 
calibration given in Table [16l If we consider the case with p = 0.5, the Monte Carlo 
method gives us the payoff expected value as 292.03 bps, with a Monte Carlo error of 
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about 1.67 bps. Again, recalling that 5*0 = 20 we can consider 292.03/20 ~ 14.6, and 
compare with CDS payoff bid ask values, as in Table [2l We see that we are within the 
7y and lOy CDS bid ask spreads but not within the ly, 3y and 5y spreads. Also in this 
case we tried to see which value for X in the ATIP model with p = 0.5 would produce 
an expected payoff close to 290. We obtained that a spread of X = 17.3 would give 
an expected payoff value of 289, so that we see that the difference between ATIP and 
SVBATIP in terms of ATIP spread is about 17.3-14.2= 3.1 bps. By comparing with 
bid-ask spreads in CDS rates R, as given in Table [T|, we see again that this difference is 
inside the ly, 7y and lOy spreads on the R quotes (which are respectively of 5, 8 an 10 
bps). Also, even if a little larger, it is comparable to the 3y and 5y spreads (both of 2 
bps), which are very narrow being related to the most traded CDS maturities. Thus we 
see that the difference is nearly negligible. 

In general we have deviated further, since we started from the payoff expected value 
under ATIP and the 165.5 expected payoff value under SBATIP. Since ATIP, SBATIP 
and SVBATIP are calibrated to the same CDS data up to five years (but SVBATIP 
is also calibrated to 7y and lOy CDS), we are seeing here that the different dynamics 
assumptions in the three models lead to different counterparty risk valuations in the 
equity return swap. The difference is not large when compared to bid ask spreads of 
CDS. We may expect more significant deviations in hedging. We have to keep in mind 
an important consideration, though. SVBATIP is calibrated not exactly and not only 
on ly, 3y and 5y CDS as the earlier ATIP and SBATIP models. Probably some of 
the difference between the price obtained with the SBATIP and the SVBATIP is to be 
attributed to this fact. 

Remark 4.2. (Nested calibration with ATIP but not with S(V)BATIP). A 

further important remark is the "nested calibration" aspect. With ATIP the calibra- 
tion is "nested", in that adding one CDS with a larger maturity does not change the 
earlier a parameters found with the calibration up to that point. In a way, this is a 
"cascade calibration" . This may be helpful with sensitivities and bucketing. Instead, in 
the S(V)BAT1P model the parameters assume a "global" role: if we add a CDS quote 
and recalibrate, all the parameters change again. In principle, every CDS quote has an 
impact on all the S(V)BAT1P model parameters at the same time. This is less desirable 
in computing sensitivities to market inputs and can lead to numerical problems, especially 
with SVBATIP where a further problem is that the calibration is not exact. 

We will address these matters in future work. 

5 Conclusions 

In general the link between default probabilities and credit spreads is best described 
by intensity models. The credit spread to be added to the risk free rate represents 
a good measure of a bond credit risk for example. Yet, intensity models present some 
drawbacks: They do not link the default event to the economy but rather to an exogenous 
jump process whose jump component remains economically unexplained. Also, when 
deahng with default correlation, a copula function must be introduced between the jump 
processes thresholds in a way that has no clear immediate relation with equity correlation, 
a source of correlation that can be used for practical purposes. 
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In this paper we introduced an analytically tractable structural first passage model 
(SVBATIP) based on scenarios on the value of the firm volatility and on the default 
barrier. This model allows for a solution to the above points. In this model the default 
has an economic cause, in that it is caused by the value of the firm hitting the default 
safety barrier value, and all quantities have a clear economic interpretation. Also, when 
dealing with multi-name products, the model allows for the introduction of the corre- 
lation in a very natural way, by simply correlating shocks in the different values of the 
firms by means of equity correlation. 

We showed how to cahbrate the model parameters to actual market data: Starting 
from CDS quotes and the scenario based barrier model, we first devised a calibration 
method based on the kernel of a matrix built by CDS with different maturities and under 
different scenarios. This approach involves a model with more parameters than market 
quotes, so that part of the parameters have to be fixed arbitrarily. After showing the 
limits of this approach, we considered a number of parameters equal to the number of 
quotes and calibrated via an optimization. This approach, examined in a case study, 
pointed out that the scenario based barrier parameterization is not describing well the 
implicit dynamics and debt structure given in CDS quotes for different maturities. An 
improvement is obtained when allowing for different scenarios also on the firm value 
volatility and also when introducing weights in the objective function, leading to a cali- 
bration error essentially within the bid ask spread. Pricing of counterparty risk in equity 
payoffs shows a partial consistency with the deterministic coefficients ATIP model in 
Brigo and Tarenghi (2004) calibrated to the same data. 

Finally, the S(V)BAT1P model can be used to build a relationship between the firm 
value V and the firm equity S (perceived as a suitable barrier payoff in terms of V itself), 
for example along the hues of Jones et al (1984) and Hull, Nelken and White (2004). 
This approach can be used to price an equity default swap. To do so we need to find 
an expression for the debt (and thus the equity) within the chosen structural model. 
Debt and equity expressions are known in closed form for time-constant and standard 
(exponential) barrier Black Cox models (see for example Bielecki and Rutkowski (2001), 
Chapter 3). Under SVBATIP we simply obtain a hnear combination of said formulas 
under each scenario and we have a closed form expression for the equity. Then we can 
price an equity default swap by means of Monte Carlo simulation of the firm value, from 
which, scenario by scenario, we deduce analytically the equity value by means of the 
found formula. This is currently under investigation. 
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